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Abstract 


Protected areas are defined as a pattern of land use in the environmental planning process. 
However, they are excluded from any physical exploitation or their use is conditional. They are 
no less important than other lands. This study aimed to evaluate the ecological degradation of 
Kosalan protected area using remote sensing and GIS techniques. Images of two time periods 
1989 and 2020 were prepared. For this purpose, after radiometric and atmospheric corrections 
of NDVI index, changes of vegetation in the region in the two time periods were evaluated. 
Then, to model the ecological degradation, eight criteria including distance from the village, 
distance from the road, landslide points, erosion intensity, slope, direction, altitude and 
vegetation were used. The criteria were weighted by hierarchical analysis process method and 
standardized by fuzzy model. Vegetation changes at different thresholds were analyzed. 
According to the vegetation situation in the region, 3 thresholds of 0.1-3, 0.1-1 and 0.1-1 were 
used. The results of this evaluation showed that the quality and density of vegetation has 
decreased a lot during 31 years. The results of ecological degradation modeling showed that the 
most effective criterion in causing degradation is distance from the village. In general, 50% of 
the area has a high potential for ecological degradation. The results showed that the greatest 
potential for ecological degradation was in the central areas of the southeast and center of the 
region, where the slope and altitude are high. The impact of slope is so high in these places. 
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این تحقیق با هدف ارزیابی تخریب اکولوژیکی منطقه حفاظت‌شده کوسالان در شهرستان مریوان در استان 
کردستان با استفاده از فن‌آوری‌های سنجش ازدور و GIS‏ انجام شده است. برای این منظور تصاویر دو 
دوره زمانی مربوط به سال‌های ۱۹۸۹ و ۲۰۲۰ تهیه گردید. پس از انجام تصحبحات رادیومتریک و 
اتمسفریک. تغییرات شاخص NDVI‏ که بیان کنندۀ تغییرات پوشش گیاهی در دو بازه زمانی موردبررسی 
است. مورد ارزیابی قرار گرفت. به‌منظور مدل‌سازی تخریب اکولوژیک از ٤‏ معیار کلی پوشش (ALS‏ 
زمین‌شناسی. فیزیوگرافی و فعالیت‌های انسانی و ۸ زیر معیار شامل فاصله از روستاء فاصله از جاده 
محدوده‌های زمین‌لغزش» شدت فرسایش. شیب. جهت. ارتفاع از سطح دریا و پوشش گیاهی استفاده 
گردید. زیر معیارها با مدل Fuzzy‏ استاندارد شدند و با روش AHP‏ وزن‌دهی گردیدند. تغییرات پوشش 
گیاهی در آستانه‌های مختلف شاخحص NDVI‏ شامل سه آستانه ۰/۳-۱ ۰/٤-۱‏ و ۰/۱-۱ مورد ارزیابی 
قرار گرفت. نتایج این ارزیابی نشان می‌دهد که کیفیت و تراکم پوشش گیاهی دارای روند کاهشی در طول 
این ۳۱ سال بوده است. نتایج مدل‌سازی پتانسیل تخریب اکولوژیکی OLE‏ داد مژثرترین معیار در ایجاد 
تخریب» فاصله از روستا است. به‌طورکلی ۵۰ درصد از مساحت منطقه دارای پتانسیل بالای تخریب 
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اکولوژیکی است. نتایج این مطالعه نشان داد استفاده از مدل‌های دقیق و تجزیه‌وتحلیل نتایج آن‌ها در محیط 
5 می تواند روش دقیقی در برآورد پتانسیل تخریب اکولوژیکی مناطق حفاظت شده باشد. 
کلیدواژه‌ها: NDVI AHP‏ فازی. مدل‌سازی» تراکم پوشش گیاهی. 


۱- مقدمه 

مناطق حفاظت شده به‌عنوان نوعی الگوی کاربری زمین در فرآیند برنامه‌ریزی محیط‌زیست تعیین شده است. 
اگرچه بهره‌برداری فیزیکی از مناطق حفاظت شده انجام نمی‌شود و یا استفاده از آن‌ها ب‌صورت مشروط صورت 
می گیرد. اهمیت آن‌ها کمتر از ارزش‌های ple‏ کاربری‌های دیگر نیست؛ زیرا این مناطق پشتوانه حمایتی مهمی برای 
بسیاری از زمینه‌ها ازجمله حفظ تنوع زیستی و توسعه پایدار کشاورزی به شمار می‌روند (روشن» ۱۳۹۷). با توجه به 
افزایش تقاضا و تشدید کاربری زمین توسط انسان‌هاء انتظار می‌رود که تغییرات کاربری زمینء مهم‌ترین عامل ایجاد 
تغییر در تنوع زیستی در آینده باشد (والز؛ ۲۰۱۱) ارزیابی آسیب‌پذیری اکولوژیکی " ناشی از اختلال مستمر انسان در 
فرآیندهای طبیعی و محیطی, توجه زیادی را در جهان به‌ویژه در زمینه تغیبرات آب‌وهوایی به خود اختصاص داده 
است (فراسر" و همکاران» ۲۰۱۱). امروزه فن‌آوری سنجش‌ازدور به‌عنوان جایگزینی برای روش‌های نمونه‌برداری 
سنتی پر هزینه و وقت‌گیر به کار می‌رود (اسلونکر* و همکاران» ۲۰۱۰؛ محمود؟ و همکاران, ۲۰۱۷؛ منزو" و 
همکاران, ۲۰۱۷). مطالعاتی در این زمینه انجام شده که می‌توان به پژوهش ادی" و همکاران (۲۰۱۷) اشاره کرد که به 
بررسی ترکیب سنجش ازدور و دانش زیست‌محیطی محلی برای نظارت بر پویایی مرانع فرفیزستان پرداختند. 
شاخحص‌های پوشش گیاهی (مانند (NDVI)‏ یا شاحص پوشش گیاهی اختلاف نرمال شده) که از تصاویر مناطق 
دور افتاده به دست آمده بودند اغلب برای نظارت بر مراتع مورداستفاده قرار گرفتند. نتایج آن‌ها نشان داد که دانش 
بوم‌شناسی محلی یک منبع با ارزش از اطلاعات و مکملی برای فن‌آوری سنجش‌ازدور است. جینگ" و همکاران 
(۲۰۲۰) با استفاده از تصاویر ماهواره‌ای روند تغییرات پوشش گیاهی منطقه با شاخص NDVI‏ را مورد ارزیابی قرار 
دادند. در این مقاله آن‌ها از تصاویر LandsatS (TM?)‏ در سال ۲۰۰۷ و تصویرگر زمین عملیاتی Landsat8‏ 


(011)در سال ۲۰۱۳ و ۲۰۱۲ استفاده کردند تا شاحص‌هایی از قبیل سبز بودن» مرطوب بودن, گرماو خشکی را 
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استخراج کنند که منعکس‌کننده کیفیت زیست‌محیطی منطقه شک بود. نتایج نشان داد که سبز بودن و رطوست 
تأثیرات مثبتی بر محیط‌زیست ia do‏ درحالی‌که گرما و حشکی اثرات منفی دارند. اکس یو و همکاران (۲۰۱۸) به 
بررسی پیش‌بینی اثرات زیست‌محیطی جمعیت بالقوه و سطح غیرقابل نفوذ با استفاده از شاخص اکولوژیکی مبتنی بر 
سنجش ازدور (RSED)‏ در ger‏ پرداختند. پیش‌بینی این تأثیر به اطلاع‌رسانی در مورد برنامه‌ریزی شهری و تصمیم- 
گیری در مورد رشد جمعیت. افزایش سطح غیرقابل نفوذ (IS)‏ و اثرات زیست‌محیطی مرتبط کمک کرد. جیانگ "و 
همکاران (X M)‏ روش ارزیابی آسیب‌پذیری زیست‌محیطی در اثر تغییر آب‌وهوا بر اساس سنجش‌ازدون در چین 
مطالعه موردی انجام داده‌اند. روش پیشنهادی مبتنی بر تعریف چارچوب IPCC*‏ و سنجش ازدور بوده است. نتایج 
نشان داد که آسیب‌پذیری زیست‌محیطی در چین متوسط یا زیاد است و تغیبرات منطقه‌ای مشخصی در توزیع 
مکانی دارد. شن * و همکاران (۲۰۱۹) مطالعه‌ای با عنوان سنجش ازدور و ارزیایی تخریب اکولوژیکی تالاب در چین 
انجام داده‌اند سنجش ازدور برای ردیابی روند تخریب اکولوژیکی در طی سال‌های ۲۰۰۰ تا ۲۰۱۵ از سه جنبه و در 
سه سطح استفاده شد. برای ارزیابی میزان تخریب تالاب از لحاظ کمی» از یک فرآیند سلسله مراتبی (AHP) |. Los‏ 
پر اساس پنج شاخص استفاده شد: نرخ تغییر منطقه تالاب › بعد JUST‏ چش‌انداز » تکه‌تکه شدن چشم‌انداز» 
زیست‌توده گیاهی و پوشش گیاهی. به‌طور حلاصه. شدیدترین تخریب تالاب در سال ۲۰۰۵ بود. در طی دوره 
۲۰۰۰-۸۵ الگوی تخریب مکانی تالاب 20188 از شمال شرقی به جنوب غربی پیشرفت کرده بود. هلوی" و 
همکاران در سال ۲۰۱۹ در مطالعه_ی از منطق فازی برای ترسیم نقشه برای میزان تولید هوموس در مناطق 
کوهستانی ایتالیا استفاده کردند. نتایج آن‌ها نشان داد که منطق فازی برای ترسیم نقشه‌های مبتنی بر تعداد داده کم و 
دارای عدم قطعیت و در مناطق ناهمگن مناسب است. تان“ و همکاران (۲۰۱۹) به بررسی شاخص ارزیابی ریسک 
زیست‌محیطی جدید و مستقیم برای تخریب محیطزیست بر اساس رویکرد منحنی پاسخ و داده‌های سنجش زدور 
پرداختند؛ و بیان کردند که اکوسیستم‌ها با تشدید فعالیت‌های انسانی و تغییرات آب‌وهوایی در جهان در معرض خطر 
تخریب محیطزیست قرار دارند. آرخی و همکاران (۱۶۰۰) به بررسی کاهش پوشش جنگلی شهرستان سردشت در 
استان آذربایجان غربی از طریق تصاویر سنجنده OLI SETM+ MSS‏ مربوط به سال‌های ۰۱۳۵۲ ۱۳۷۹ و ۱۳۹۷ 
پرداختند. تصاویر ذکر شده مورد پیش‌پردازش و پردازش قرار دادند و به دو طبقه جنگل و غیرجنگل طبقه‌بندی 
کردند. نتایج نشان داد که در طول ۶۱ سال حدود ۳۳۷۲۱ هکتار از سطح جنگل‌های شهرستان سردشت کاسته شده 
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است. با توجه به نتایج مشخص شد که از متغیرهای توپوگرافی» پارامترهای فاصله از جاده و فاصله از روستا بر ميزان 
تخریب جنگل بیشترین تأثیر را داشته‌اند. همان‌طور که مطالعات انجام شده OUS‏ داد فن‌آوری سنجش ازدور و بررسی 
معیارهای مختلف و مدل‌سازی آن به کمک مدل‌های چند معیاره مانند AHP‏ و منطق فازی می‌تواند در شناسایی 
تخریب اکولوژیک مناطق مختلف مؤثر باشد. لذا در این تحقیق مدل‌سازی تخریب اکولوژیک در محیط سامانه‌های 
اطلاعات جغرافیایی با استفاده از مدل‌های منطق فازی و AHP‏ انجام برای dikes‏ حفاظت شده کوسالان انجام شد. 


-Y‏ مواد و روش 

۱-۲- منطقه مورد مطالعه 

منطفه مطالعاتی در غرب ایران» استان کردستان و در شهرستان مریوان قرار گرفته است. کوسللان در محدوده 
طول‌های جغرافیایی ۲۷ ۱۱۳ 47 تا ۳۰۳۲۷ EV‏ شرقی و عرض‌های جغرافیایی Yo VV"‏ تا Yo YY OW‏ 
مریوان (مابین شهرهای سروآباد- مریوان و پاوه) قرار گرفته است (شکل A)‏ هم‌چنین بخشی از این منطقه در استان 
کرمانشاه واقع شده است که این مطالعه بخش کوسالان را در بر گرفته که مساحت آن TV‏ هزار هکتار است در این 
منطقه تاکنون ۱۱۷ گونه پرنده» ۲۳ گونه پستاندان ۱۷ گونه خزنده و ۲۶۲ گونه گیاهی شناسایی شده است. مهم‌ترین 
پوشش گیاهی این منطقه شامل جنگل‌های بلوط است. گونه‌های شاخص جانوری شامل: سیاه‌گوش کل و بز 
خرس قهوه‌ای. پلنگ» کرکس, هما و عقاب طلائی است (دیده بان محبط‌زیست و حیات وحش. ATAN‏ 
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East Azarbalján 
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ارزبابی تخربب اکوئوژبکی منطقه حفاظت شسده 
کوسالان با استفاده از سنجش از دور و GIS‏ 


-Y-Y‏ داده‌ها 


انواع داده‌های مورداستفاده در این مطالعه و jsi‏ و مشخصات آن‌ها در جدول شماره ۱ نشان داده شده است. در 
این مطالعه نقشه روستاهاء حاده‌ها؛ تقاط زمین‌لغزش» فرسایش و DEM!‏ از سازمان نقشه‌برداری کشور ag‏ گردید» 
همچنین EET‏ ماهواره‌ای منطقه از وب‌سایت سازمان زمین‌شناسی LS yl‏ دانلود و تصحیحات مربوطه — 
استخراج اطلاعات بر روی آن‌ها انجام شد (جدول ۱). 


حدول ۱- داده‌های مورداستفاده و مشخصات آن‌ها 


داده مأخذ داده فرمت داده مقیاس داده دوره زمانی سال اخذ داده 

۱۳۹ سازمان نقشه‌برداری رس‌تر ۰ متر‎ DEM 

تصاویر وب‌سایت زمین‌شناسی (SRTMDG, ol‏ رس‌تر ۰ متر ۱۹۸۹-۶۰ ME‏ 
پوشش گیاهی سازمان جنگل‌ها و مراتع برداری (وکتور) ۱/۵۰ ۱۳۹۹ 


1 Digital Elevation Map 
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داده مأخذ داده فرمت داده مقیاس داده دوره زمانی سال اخذ داده 
روستاها سازمان نقشه‌برداری برداری (وکتور) ۱۵۰۰۰ ۱۳۹۹ 
جاده‌ها سازمان نقشه‌برداری برداری (وکتور) ATTEND‏ ۱۳۹۹ 
نقاط زمین‌لغزش سازمان نقشه‌برداری برداری (وکتور) ۱۰۰۰ ۱۳۹۹ 
فرسایش سازمان نقشه‌برداری برداری (وکتور) ۱۵۰۰۰ ۱۳۹۹ 


as e‏ امین 

۱-۳-۲- پایش تغییرات پوشش گیاهی 

برای ارزیابی تغییرات پوشش گیاهی از شاحص | NDVI‏ استفاده شد. ابتدا تصاویر دو دوره زمانی ۱۹۸۹ و ۲۰۲۰ 
از ماهواره منابع زمینی Landsat‏ دانلود شده و پس از انجام تصحیحات رادیومتریک و اتمسفری با استفاده از روش - 
های (FLAASH -Drak subtract)‏ با شاعص NDVI‏ تغییرات پوشش گیاهی دو دوره بررسی شد. به‌منظور 
کنترل کیفیت داده‌های مورداستفاده و LAIST‏ از اینکه خطاهای سیستماتیک و غير سیستماتیک تا چه حد در تصحیح 
سین bos‏ لی با At‏ اس Tani‏ نی 5 رس coal,‏ قرار SS‏ مس از تما اناه ورت 
تک باندی و ترکیبات رنگی مختلف با استفاده از نرم‌افزارهایی که قابلیت نمایش و بزرگنمایی بالایی دارند. به کمک 
بزرگنمایی آن‌ها در بخش‌های مختلف تصاویر از لحاظ خطاهای رادیومتری نظیر راه راه شدگی» حطای زیر هم 
قرارگرفتن خطوط اسکن. پیکسل‌های دوبله و همچنین اثرات اتمسفری نظیر ابر موردبررسی قرار گرفتند (نظم فر و 
همکاران ۱۶۰۰). 


حدول۲- مشخصات تصاویر ماهواره‌ای مورداستفاده 


دوره زمانی ماهواره سنحنده ماه رزولوشن 
TM Landsat ۱۹۸۹‏ اواخر ژوئن ۰ متر 
قاط OLI Landsat‏ اواخر ژوئن ۳۰ متر 


تصاویر دو دوره زمانی با فاصله ۳۱ سال در ماه ژوئن (کیفیت تصویر در این ماه به علت وجود پوشش PLS‏ و 


همچنین عدم وجود ابر بهتر است) دانلود و تصاویر هر دو دوره در یک ماه برداشت شد تا تغییرات پوشش گیاهی 
در تراکم مشابه در دو دوره مختلف. به‌خوبی قابل قیاس باشند (جدول AY‏ 
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۷-۳-۲- مدل‌سازی تخریب اکولوژیکی 


در این بخش با استفاده از ۸ زیر معیار به ارزیابی تخریب اکولوژیکی پرداخته شد. ابتدا زیر معیارها در محیط 
نرم‌افزار 10.7 ArcGIS‏ ماده‌سازی شدند و سپس با استفاده از منطق فازی استاندارد شدند. برای وزن‌دهی به معیارها 
3.3 فطل سیب Siw‏ تفای ان رازفا یرای طق رن oa‏ این سای شاه فان اين 
نیز در سایر تحقیقات مورداستفاده قرار گرفته بودند که با بررسی وضعیت جغرافیایی منطفه آن‌ها حذف شدند. 
معیارهای مهمی که از مدل‌سازی حذف شدند شامل اقلیم و تغییر کاربری اراضی بودند که معیار اقلیم با توجه به 
اقلیم OLS‏ حذف گردید. تغییر کاربری نیز با استفاده از طبقه‌بندی نظارت شده تصاویر ماهواره لندست و همچنین 
استفاده از GL.‏ گوگل ارث در دامن زمانی مطالعة حاضر موردبررسی قرار گرفت که در هر دو روش به دلیل 
کوهستانی بودن و توپوگرافی بسیار سخت Of‏ تغییر کاربری اراضی مشاهده نشد. سایر معیارها با توجه به اهمیت 


آن‌ها و کاربرد در مدل تخریب مورداستفاده قرار گرفتند. 


حدول ۳- معیارهای مورداستفاده در مدل‌سازی تخریب اکولوژیک 


ردیف معیار زیرمعیار منبع 
۱ پوشش پوشش گیاهی NDVI‏ سپهر و همکاران» ۱۳۹۶ 
۱ فاصله از bE‏ زمین‌لغزش شیرمحمدی و همکاران ۱۳۹۵ 
iie i‏ فرسایش سپهر و همکاران ۱۳۹۶ 
شیب سپهر و همکاران ۱۳۹۶-شیرمحمدی و همکاران ۱۳۹۵ 
۳ فیزیوگرافی جهت سپهر و (OU‏ ۱۳۹۶- شیرمحمدی و همکاران.» ۱۳۹۵ 
ارتفاع سپهر و همکاران. ۱۳۹۶- شیرمحمدی و همکاران ۱۳۹۵ 
فاصله از روستا شیرمحمدی و همکاران ۱۳۹۵ 
3 فعالیت انسانی 
فاصله از جاده شیرمحمدی و همکاران ۱۳۹۵ 


۱-۲-۳-۲- فازی‌سازی لایه‌ها در سیستم اطلاعات جغرافیایی 

هدف از استاندارسازی Lag‏ همسان نمودن آن‌ها برای به‌کارگیری در سامانه پشتیبان تصمیم‌گیری مورداستفاده 
است که در این تحقیق مدل تخریب اکولوژیک است. در این مدل‌سازی معیارهای مورداستفاده با منطق فازی 
استانداردسازی شد. با استفاده از مبانی تئوری فازی» که هر موضوع یا مفهومی دارای یک درجۀ عضویت است؛ 


FA‏ جغرافیا و مخاطرات محیطی شمارة اول 


جدول E‏ لایه‌ها و alg‏ فازی مورداستفاده برای استانداردسازی آن‌ها و هم‌چنین نقاط کنترلی (control point)‏ و 


نقطه وسط (midpoint)‏ نشان داده شده است. 


حدول ۶- استانداردسازی لایه‌های مدل تخریب اکولوژیکی 


ردیف معیار زیرمعیار تابع فازی نقطه کنترلی ۱ نقطه وسط نقطه کنترلی ۲ 
۱ پوشش پوشش گیاهی NDVI‏ خطی افزایشی y KA‏ 
Y‏ فاصله از نقاط زمین‌لغزش حطی کاهشی ١کم‏ 
T My‏ فرسایش خطی افزایشی ۱ ۳ 
t‏ شیب خطی افزایشی ۱ 8 
۵ فیزیو گرافی جهت Near‏ 10 
1 ارتفاع خطی افزایشی ۸.۰ Ye‏ 
۷ - فاصله از روستا خطی کاهشی ٦‏ کم 
u ^‏ فاصله از جاده خطی کاهشی ١کم‏ 


برای تعدادی از معیارها از تابع عضویت فازی خطی استفاده شده است. این معیارها با توجه به ماهیت لایه 
جغرافیایی و هم‌چنین ارتباط با مدل اکولوژیکی با این تابع استاندارد شده‌اند. به‌عنوان‌مثال gY‏ فاصله از روستاء از تابع 
خحطی کاهنده استفاده شده است. بدین معنی که هر چه فاصله از روستا بیشتر باشد ميزان پتانسیل تخریب اکولوژیک 
کمتر است و هر چه فاصله کمتر باشد پتانسیل تخریب بیشتر است؛ که این با eU‏ خطی کاهنده مطابقت دارد. نقاط 
کنترلی نیز با توجه به تابع فاصله اقلیدسی به دست آمده است که از صفر شروع می‌شود و تا دورترین نقطه مرز 
منطقه از روستاهای موجود است. بیشترین فاصله از روستا تا مرز منطقه کوسالان. شش کیلومتر است. این تابع برای 
هفت زير معیار به‌صورت کاهش و افزایشی استفاده شده است. برای زیر معیار جهت شیب از تابع عضویت نزدیک 
(NEAR)‏ استفاده شده است که تابعی شبیه تابع گوسین یا توزیع نرمال است و دارای نقطه وسط است. با معرفی 
نقطه وسط به‌عنوان با ارزش‌ترین نقطه یا با بیشترین پتانسیل تخریب در این مطالعه. هر چه از این نقطه به طرفین 
فاصله بگیریم ارزش يا پتانسیل تخریب OF‏ کاسته می‌شود. با توجه به اینکه جهت جغرافیایی با ميزان دريافت تابش 
خورشید و میزان تبخیر و سپس آب قابل دسترس برای پوشش گیاهی و رطوبت ارتباط مستقیم داردء بنابراین نقطه 
۵ درجه که جنوب غربی هست به‌عنوان بالاترین پتانسیل تخریب به تابع نزدیک معرفی شده است و نقطه وسط 
۱ نیز به‌عنوان میزان پخش‌شدگی Of Spread)‏ که بیشترین میزان پخش‌شدگی در این تابع است. معرفی شد. 


توابع عضویت تمامی معیارها در جدول ۶ آمده است. 
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در شکل شماره n Y‏ خروجی mE‏ که در نرم‌افزار GIS‏ استانداردسازی شده با مدل فازی AE adal‏ 
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شکل ۲- aca‏ فازی معیارهای استفاده شده 


۲-۳-۲-۲- وزن‌دهی به معیارها با روش فرایند تحلیل سلسله مراتبی AHP‏ 


جهت بیان اهمیت هر معیار نسبت به معیارهای دیگر Lb‏ این معیارها وزن‌دهی شوند. در ادبیات تصمیم چند 


معیاری» روش‌هایی نظیر رتبه‌بندی, درجه‌بندی مقایسه دو به دو و تحلیل موازنه- جایگشتی در وزن‌دهی معیار بر 


al‏ قضاوت‌های تصمیم گیران ارائه شده است. درصورتی که زمان و هزینه کم و به دست آوردن وزن‌هابه طریق 
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آسان مهم‌ترین مسئله باشد» یکی از روش‌های رتبه‌بندی یا نسبتی توصیه می‌شود؛ اما هنگامی که دقت. اساس کار 
مسئله باشد» روش مقایسه دوتایی و با محبط زیست " تحلیا توازن مناسب است L)‏ ^ لجفسک, ù‏ ۵ در اسن مطالعه 
پرسشنامه‌ها بین متخصصان منابع طبیعی و محیطزیست پخش شد و نتایج با نرم‌افزار Export choise‏ مورد 


تجزیه‌وتحلیل قرار گرفت که در جدول شماره ۵ نشان داده شده است و در OLL‏ نرخ ناسازگاری VAR‏ به دست آمد. 


جدول ۵- اوزان معیارها 


ردیف نام معیار وزن بر اساس روش AHP‏ 
۱ فاصله از روستا 2 
۲ پوشش گیاهی UYNW‏ 
۳ فاصله از جاده ۰۱۹ 
3 نقاط زمین‌لغزش TIN‏ 
۵ فرسایش ۰/4 
1 شنیب MA‏ 
۷ جهت ۰۳ 
۸ ارتفاع 4 


در مرحله آخر برای تولید نقشه نهایی پتانسیل تخریب لایه‌های تولید شده در نرم‌افزار سامانه اطلاعات جغرافیایی 
روی هم‌گذاری می‌شوند. 
-Y‏ نتایج و co‏ 


همان‌طور که در شکل ۶ مشاهده می‌شود. میزان پوشش گیاهی برای دو آستانه ۰/۲-۱ و ۰/1-۱ در سال ۱۹۸۹ 
نشان داده شده است. آستانه ۰/٤-١‏ پوشش با کیفیت را نشان می‌دهد و آستانه ۰/۲-۱ پوشش با کیفیت کمتر از 


آستانه ۱/۶-۱ را نشان می‌دهد (براتکومار ' و محمد (Yo‏ 


( 
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شکل ۳- پوشش گیاهی (شاخص (NDVI‏ در سال ۱۹۸۹ 


ارربابی تریب اکولوزیکی منطقه حفانظت ده 
کوسالان با استفاده از سسنجتش از هور و CIS‏ 


tant, 


-Y col »NDVI 8 
ras o NDVI RES 
ar CO 


NDVI شاخص‎ 
0.4 +0.3 a di بر اسلی‎ 
2020 
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شکل £- پوشش گیاهی (شاخص (NDVI‏ در سال ۲۰۲۰ 
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ارزبابی تخرهب اکولوزیکی in dica‏ شسده 
کوسالان ما فده از سنجنی از هور و 615 


[n 
0.1-1 + NDVI 1989 7 
0.1-1 ت‎ NDVI 2020 0 


KOSALAN PROTECTED AREA 
UTM ZONE 38 N 


x 


در شکل شماره۳ و 6 نقشه پوشش گیاهی بر اساس دو آستانه ۰/۲-۱ و ۰/٤-۱‏ در شاخحص NDVI‏ به ترتیب 
در سال‌های ۱۹۸۹ و ۲۰۲۰ نشان oala‏ شده است. در شکل شماره ۵ نتیجه شاخص NDVI‏ در دامنه ۰/۱-۱ در هر 
دو دوره را OUS‏ می‌دهد. همان‌طور که مشاهده می‌شود جز قسمت‌های کوچکی که شامل اراضی ple‏ و روستاها 
می‌شود بقیه مناطق دارای پوشش می‌باشند. پوشش هر دو دوره بر هم منطق شده‌اند و در هر دو دوره در این آستانه 
کل منطقه بر اساس شاخص NDVI‏ دارای پوشش می‌باشند. نتایج پوشش گیاهی با شاخص NDVI‏ نشان داد که 
در طول زمان موردبررسی کیفیت پوشش گیاهی تغییرات زبادی داشته است. با استفاده از این شاخص می‌توان 
وضعیت اکولوژیکی منطقه را از لحاظ پوشش ALS‏ بررسی نمود. در صورت کاهش سطح پوشش ALS‏ و جنگلی 
منطقه می‌توان پی برد که وضعیت اکولوژیکی منطقه نیز در حال افت است. استفاده از شاخص NDVI‏ برای نشان 
دادن تغییرات اکولوژیکی در محط‌های مختلف توسط محققین زیادی بررسی شده است (ژانگ" و همکاران, ۲۰۲۲؛ 
پراوالی " و همکاران ۲۰۲۲؛ لی" و همکاران, ۲۰۲۱؛ جیانگ" و همکاران» ۲۰۲۱). نتایج این مرحله از مطالعة حاضر با 
نتایج گیلسپای"* و همکاران در سال ۲۰۱۸ همخوانی دارد گیلسپای و همکاران در سال ۲۰۱۸ نیز در مطالعۂ از فن- 
آوری سنجش‌ازدور برای OLS‏ دادن تغییرات NDVI‏ در Ais‏ حفاظت شده در جنوب کالیفرنیا استفاده کردند. نتایج 


1 Zhang 

2 Pravalie 
3 Li 

4 Jiang 

5 Gillespie 


oF‏ جغرافیا و مخاطرات محیطی شمارة اول 


آن‌ها نیز کاهش NDVI‏ را در سطح منطقه OUS‏ داد و همچنین آن‌ها بیان داشتند که فن‌آوری سنجش‌ازدور برای 
بررسی تغیبرات پوشش گیاهی قابل‌اعتماد است. ازآنجایی که محدوده طیفی Landsat 8 OLI‏ باریک‌تر از Landsat‏ 
TM‏ 4-5 است. NDVI‏ مشتق شده از این ماهواره‌ها مقادیر متفاوتی دارد. باریک بودن محدوده طیفی سنجنده OLI‏ 
سبب می‌شود که مقدار NDVI‏ را در مقایسه با ماهواره Landsat 4-5 TM‏ بیشتر OUS.‏ دهد. با توجه به این امر که 
یکی از تصاویر مورداستفاده در این مطالعه تصویر Landsat 4-5 TM‏ است (سال OM‏ و تصویر دیگر تصویر 
ماهواره JU) Landsat 8 OLI‏ ۲۰۲۰) است؛ لذا برای حل این مشکل» شاخحص‌های این دو سنسور بايد e JUS‏ شوند 
(Gooshbor et al., 2016)‏ اما با توجه به نتایج به دست آمده از شاحص NDVI‏ در سال ۲۰۲۰ که از نظر دامنه 
عددی کمتر از سال ۱۹۸۹ است می توان نتبجه گرفت که کاهش کیفیت اکولوژیک منطقه شدید بوده و با توجه به 


عدم افزایش شاحص NDVI‏ در سال ۲۰۲۰ نیازی به انجام کالیبراسیون نیست. 


۲-۳- پهنه‌بندی پتانسیل تخریب اکولوژیک 

در این بخش لایه‌های استاندارد شده مرحله قبلی همراه با وزن آن‌ها روی‌هم گذاری شده و نقشه پتانسیل 
تخریب اکولوژیکی منطقه حاصل شد. این کار با دستور map algebra‏ و با کد زیر اجرا می‌شود. 

+ GIA +" فازی شیب‎ Y) + )۰/۰۸۹ S ("لایه فازی ارتفاع‎ + G/VW فازی پوشش گیاهی"*‎ a") 
+ )۰/۰۹۹ ag" فازی فرسایش‎ YU) + )۰/۰۹۳ P ("لایه فازی جهت شیب‎ + C/V OV * فازی فاصله از جاده‎ Y") 


)۰/۱۸۹ فازی فاصله از روستا"٭‎ a") + )۰/۱۰۹ *" فازی نقاط زمین‌لغزش‎ NU) 


از نمی بویت اکولوزبکی متخت PM‏ 
کومالان استفاده ار مهن از دور و OS‏ 


Por sae | 
پتانسیل تخریپ اکولوژیکی‎ 
eu 


ip = 
IAOW ابی‎ 
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wre ۳ 


شکل *- نقشه پتانسیل تخریب اکولوژیکی 
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در این مطالعه با توجه به روش استانداردسازی بر اساس ماهیت تحقیق, هر چه ارزش سلول SVL‏ باشد پتانسیل 
تخریب بالاتر و هر چه ارزش سلول پایین‌تر باشد پتانسیل تخریب پایین‌تر است. در شکل شماره 1 نقشه پتانسیل 
تخریب منطقه نشان داده شده است. بر اساس عملکرد تابع ارزش بین ۰/۲۸ تا ۰/۷۸ به دست آمده است. مناطق با 
پتانسیل تخریب کم سمت ارزش کم بوده و مناطق با پتانسیل تخریب بالاء سمت ارزش زیاد می‌باشند. همان‌طور که 
در VS ste‏ نشان داده شده است. پتانسیل تخریب در منطقه به ۵ طبقه تقسیم شده است: 


جدول —V‏ نقشه پتانسیل تخریب اکولوژیکی 


پتانسیل تخریب اکولوژیک مساحت (هکتار) درصد 
خیلی کم jm M‏ 
کم 040 ۳/07 

Yow ۸۱ متوسط‎ 

IY Mew زیاد‎ 

vM IT خحیلی زیاد‎ 


WA‏ هکتار دارای پتانسیل تخریب اکولوژیک خیلی کې ۳۵۹۵ هکتار تخریب کم ۹۶۸۱ هکتار تخریب متوسط 
۷ هکتار تخریب زیاد و VW‏ هکتار دارای پتانسیل تخریب خیلی زیاد می‌باشند. هم‌چنین میزان پتانسیل 
تخریب به‌صورت درصدی» ۰/1۱ درصد تخریب خیلی کم ۱۳/۵۱ درصد تخریب کم» ۳۵/۷۷ درصد تخریب 
متوسط 1۳/۸ درصد تخریب زیاد و ۱۱/۰۹ درصد تخریب خیلی زیاد هست. در شکل V‏ مناطق با بالاترین پتانسیل 
تخریب به رنگ قرمز نشان داده شده است. این مناطق بیشتر در سمت شمال غرب منطقه قرار گرفته‌اند. هم‌چنین 
بخشی از این مناطق در جنوب شرق و هم‌جوار با روستاها و جاده‌های ارتباطی آن‌ها واقع شده است. مناطق سبز 
رنگ دارای کمترین و مناطق قرمز رنگ دارای بالاترین پتانسیل تخریب می‌باشند. 
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شکل ۸- نقشه پهنه‌های با پتانسیل تخریب اکولوژیکی بالا 


Juss با‎ gangga هلف از این کار‎ cathe واده‎ GL Dee Sy ye نالا‎ Se bab yak IRS در‎ 
dex yh Siete odiosa, Si Sia aet o Labeo obs aes 
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مختلف شامل مناطق بیابانی (لیو" و همکاران» ۲۰۲۲ مناطق حفاظت شده (نعمت‌اللهی" و همکاران ۲۰۲۲ لی" و 
همکاران. (THY‏ محیط‌های شهری SS)‏ و همکاران» ۲۰۲۱) و ب‌صورت کلی‌تر در کشورهای توسعه یافته 6-7 
(عثمان" و همکاران ۲۰۲۱) بررسی کرده‌اند. در همه این مطالعات کیفیت اکولوژیک کاهش داشته است و یا در 
بهترین حالت ممکن در قسمت‌هایی که دور از دسترس فعالیت‌های انسانی بوده ثابت و بدون تغییر بوده است 
به‌صورت ویژه مطالعة نعمت‌اللهی و همکاران (۲۰۲۲) نشان داد که فعالیت‌های انسانی مانند حطوط حمل‌ونقل جادۀ 
سبب شده که مناطق مناسب برای زندگی حیات‌وحش در مناطق حفاظت شده استان چهار محال بختیاری در ایران 
دارای ریسک بالا و متوسط باشد که این نتایج با نتایج پژوهش حاضر در اثرگذاری بالای فاصله از جاده در تخریب 
اکولوژیک منطقه همخوانی دارد در این راستا امامی و شهریاری (۱۳۹۸) نیز نشان دادند که عامل نزدیکی به جاده 
بیش‌ترین سهم را در ایجاد وقوع آتش‌سوزی در منطقه حفاظت‌شده ارسباران دارد. یکی دیگر از نتایج پژوهش حاضر 
اثرگذاری نزدیکی به مراکز جمعیتی یا همان روستاهای منطقه در افزایش ميزان تخریب اکولوژیکی است که این نتایج 
با نتایج گرانادوس" و همکاران در سال ۲۰۲۲ و شیر محمدی و همکاران در سال ۱۳۹۵ مطابقت دارد. گرانادوس و 
همکاران ۲۰۲۲ در مطالعة تغییرات پوشش گیاهان در معرض خطر انقراض را با استفاده از شاخحص NDVI‏ در 
مناطق حفاظت شده جنوب اروپا موردبررسی قرار دادند. نتایج آن‌ها Sly‏ عدم‌تغییر پوشش در مناطق بدون دخالت 
انسان بود. شیرمحمدی و همکاران (۱۳۹۵) در مطالعة اثرات توسعه را بر منطقه حفاظت شده کرکس بررسی کردند. 
آن‌ها به این نتیجه رسیدند که فعالیت‌های انسانی از مهم‌ترین معضلات منطقه بوده و بیشترین تأثیر را در تخریب این 
مناطق داشته است. با توجه به نتایج حاصل از این پژوهش می‌توان گفت استفاده از فناوری سنجش‌ازدور در JUS‏ 
استفاده از مدل‌های AHP Fuzzy‏ در محیط GIS‏ روش دقیقی برای برآورد پتانسیل تخریب اکولوژیکی در مناطق حفاظت 


شده است. 


۶- جمع‌بندی 

این تحقیق با هدف ارزیابی تخریب اکولوژیکی منطقه حفاظت‌شده کوسالان با استفاده از سنجش زدور و GIS‏ 
انجام شده است. ابتدا به‌منظور تغییرات پوشش گیاهی از فن‌آوری سنجش‌ازدور استفاده شد. تصاویر دو دوره زمانی 
مربوط به سال‌های ۱۹۸۹ و ۲۰۲۰ angi‏ و با استفاده از نرم‌افزار 5.3 Envi‏ مورد ارزیابی قرار گرفت. سپس به‌منظور 
مدل‌سازی تخریب اکولوژیک از ۸ زیر معیار استفاده شد. با روش فرآیند تحلیل سلسله مراتبی وزن‌دهی شدند و با 


1 Liu 

2 Nematollahi 
3 Li 

4 Feng 

5 Usman 

6 Granados 
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منطق فازی استاندارد شدند. با روی هم گذاری لایه‌ها نقشه پتانسیل تخریب اکولوژی € و موردبررسی و ارزیابی 
قرار گرفت. بر اساس شاخحص NDVI‏ (شاحص نرمال‌شده تفاوت پوشش گیاهی) که یک شاخحص SIS‏ ساده 
است که در تحلیل‌ها و اندازه‌گیری‌های سنجش‌ازدور و ارزیابی وجود یا عدم وجود پوشش lS‏ یک منطقه کاربرد 
oul‏ تغییرات پوشش گیاهی در آستانه‌های مختلف بررسی و تحلیل شد. با توجه به وضعیت پوشش ALS‏ در 
منطقه از ۳ آستانه برای این امر استفاده شد. آستانه ۰/۳-۱ ۰/٤-۱‏ و ۰/۱-۱ استفاده شد. نتایج آستانه ۱۱۰/۱ 
NDVI | az Us‏ نشان داد که منطقه موردمطالعه در طول ۳۱ سال میزان تخریب و يا تبدیل مراتع به دیم‌زار خیلی 
کمی داشته است چون تمامی منطقه در هر دو تصویر بازه‌های زمانی ذکر شده دارای پوشش ALS‏ در تراکم ۰/۱ تا 
| بود ولی نتایج آستانه‌های دیگر شامل ۰/۳-۱ و ۰/٤-١‏ که آستانه‌های بیان کننده کیفیت پوشش گیاهی بر اساس 
تراکم پوشش است نشان داد که تراکم پوشش گیاهی در طول ۳۱ سال کاهش داشته است. به نظر می‌رسد کاهش 
تراکم پوشش گیاهی منطقه به دلیل چرای مفرط و بیش از ظرفیت مراتع باشد که افزایش جمعیت انسانی روستاهای 
حاشیه منطقه و به طبع آن افزایش نیاز به استفاده از مراتع این مناطق جهت تعلیف دام سبب این امر شده است. البته 
در کنار چرای مفرط عامل تغییر اقلیم و کاهش بارش برفی منطقه نیز می‌تواند تأثیرگذار باشد که نیاز به بررسی بیشتر 
در مطالعات آتی دارد. به‌صورت IS‏ می‌توان گفت شرایط فیزیو گرافی شامل کوهستانی بودن منطقه همراه با دامن 
ارتفاعی و درجهٌ شیب بالا سبب شده که عامل مهمی مانند تغییر کاربری اراضی که سبب از بین رفتن مناطق تحت 
استفاده حیات‌وحش است در این منطقه وجود نداشته باشد. از طرف دیگر عوامل انسان‌ساخت دیگری مانند 
سکونتگاه‌های روستای و احداث جاده در مناطق iilo‏ منطقة حفاظت شده سبب افزایش پتانسیل تخریب 
اکولوژیکی این مناطق شده است. درجه شیب بالا همراه با دامنة ارتفاعی زیاد از مهم‌ترین عوامل افزایش پتانسیل 
تخریب اکولوژیک منطقه است که سبب ایجاد اکوسیستم‌های بسیار حساس و شکننده در مقابل فشارهای بیرونی 
شده است. این امر نشان می‌دهد در کنار این محدودیت‌های طبیعی هرگونه دخالت انسانی حتی به‌صورت محدود 
مانند چرای بیش از ظرفیت مرتع و یا ایجاد جاده به‌منظور توسعهٌ صنعت توریسم می‌تواند خسارات جبران‌ناپذیری 
در منطقه به بار آورد. 
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